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' Although the physics of normal metals (N) in close proximity to superconductors (S) has been studied extensively 
^ for more than thirty years, it is only in the past decade that experiments have been able to probe directly the region 
close to the NS interface at temperatures far below the transition temperature Tc of the superconductor. These 
experiments have been made possible by the availability of microlithography techniques that enable the fabrication 
• of heterostructure devices with submicron scale features. This size scale is comparable to the relevant physical length 
\ scales of the problem, and consequently, a number of new effects have been observed. A variety of systems have 
been studied, the variation primarily being in the type of normal metal coupled to the superconductor. Canonical 
' normal metals such as copper and gold, semiconductors, insulators, and ferromagnets have been employed. Although 
' ^ a variety of theoretical techniques have been used to describe proximity-coupled systems, the quasiclassical theory of 
' superconductivityiiSiii^iSiL^iSiiSiiiiiS has proved to be a remarkably powerful tool in understanding the microscopic 
, basis for the remarkable effects observed in these systems. 
Q ■ A number of excellent recent articlesi^ii^ii^iiSiiLiSiiS have explored the application of the quasiclassical theory of 
Cj ', superconductivity to proximity coupled systems. In this review, a self-contained development of the quasiclassical 
. ■ theory is presented, starting from non-equilibrium Keldysh Green's functions for normal metal systems. If the normal 
"j^ I metal is clean, quasiparticles in the normal metal travel ballistically over long length scales; in the samples studied 
in the majority of recent experiments, however, the quasiparticles are scattered elastically within a short distance, so 
that the quasiparticle motion is diffusive. Here, we shall concentrate on this diffusive case, where the elastic scattering 
, length £ is the shortest relevant length scale in the problem. We shall also restrict ourselves to the case where the 
^ ■ superconductor is of the canonical s-wave type, avoiding complications with non-spatially symmetric order parameters 
O \ within the superconductor. 

CN \ 1. TRANSPORT EQUATIONS IN THE DIFFUSIVE APPROXIMATION 
> 

. As an illustration of some of the issues that arise in dealing with non-equilibrium transport in mesoscopic diffusive 

ly-j ' systems, we consider first the classical Boltzmann equation in the diffusion approximation. Consider a one-dimensional 

fvq diffusive normal metal wire of length L. We assume that all dimensions of the wire are larger than £, the elastic mean 

T-H free path of the conduction electrons. f{E,x), the distribution function of electrons at energy E and at a point x 

, along the wire, obeys the diffusion equation 

i,!fZg:£>+C(/).o. (1.1) 
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I ' Here, D = (l/3)vi?f is the three-dimensional electron diffusion constant, with vp being the Fermi velocity of the 



electrons. C(/) is the collision integral, which takes into account inelastic scattering of the electrons, and itself 
depends on the distribution function /. If we consider a mesoscopic wire whose length L is much shorter than any 
^ inelastic scattering length Lin, this term can be set to zero. 

H^* This diffusion equation has a simple solution in some specific geometries. Consider the case of a wire of length L 

sandwiched between two normal 'reservoirs,' which we shall call the left (at x — Q) and right (at x = L) reservoirs. 
The reservoirs are defined as having an equilibrium distribution function /l(£') and fniE) respectively. Then the 
?H solution of the diffusion equation, Ean. ljl.lll . under the condition L <^ (so that C(/) = 0), and subject to the 
boundary conditions that f{E,x) ~ fL{E) at the left reservoir, and f{E,x) = fniE) at the right reservoir, is given 

by2a 

f{E,x) = [fR{E) - h{E)] I + h{E) (1.2) 
The electrical current / through the wire in the diffusion approximation is given by 

I = -eAD J N{E)^dE (1.3) 



FIG. 1: Nonequilibrium electron distribution function / in a normal wire with a potential V applied across it, as a function of 
energy E and position x. The position x is normalized to the length L of the wire. 



where A is the cross-sectional area of the wire. With f{E,x) given by Ean. H1.2|l . this can be written 

eAD 
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N{E)[fR{E)- fL{E)]dE 



(1.4) 



In order to obtain a finite current, we must apply a voltage across the wire. Let us apply a voltage V to the left 
reservoir, keeping the right reservoir at ground (1/ = 0). This has the effect of shifting the energies of the electrons 
in the left reservoir by —eV, so that the electron distribution function there is given by fhiE) = fo{E — eV), where 
fo{E) is the usual equilibrium Fermi distribution function 
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(1.5) 



and the energy E is measured from the Fermi energy Ep. Figure 1 shows the distribution function f{E, x) in the wire. 
It has a step-function form, which varies linearly along the length of the wire. Evidence for such a non-equilibrium 
distribution has recently been seen experimentally in a series of beautiful experiments by the Saclay groupn^ 
In the limit of small voltages V, the difference of distribution functions in Ean. l|1.4|) can be expanded as 



fBiE) ~ UiE) = foiE) - foiE ~ eV) « /o(i?) 



ME) + eV 
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dE 



= -eV 



dfo 
' dE 



so that 



/ = 



dE 



N{0)e^AD 



V. 



^ (1-6) 



(1.7) 



where iV(0) is the density of states at the Fermi energy. In performing the integral, we have assumed that we are at 
low enough temperatures so that the derivative of the Fermi function {—dfo/dE) can be approximated by S{E). The 
electrical conductance is then given in the Nernst-Einstein form G — N{0)e'^D{A/L) — aQ{A/L). 
The thermal current through the wire is given by an expression similar to Ean. (|1.4|) : 



AD 
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EN{E)[fu{E)- fL{E)]dE 



(1.8) 



The critical difference between the expression for the thermal current and the expression for the electrical current 
is the presence of an additional factor of E in the integrand in Ean. (|1.8|l . which has important consequences in the 
calculation of the thermal properties. 
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To obtain expressions for the thermal coefRcents, let the temperature of the left reservoir be T and its voltage 
V , and let the temperature of the right reservoir T + AT and its corresponding voltage = 0. The distribution 
function in the left electrode is then Jl (E) — fo{E — eV, T) , and the distribution function in the right electrode is 
fii{E) = fo{E, T + AT). In the hmit of small AT, V, the difference in the distribution functions in Ean. l|1.4f) can be 
expanded as 



fR{E)-fL{E) = 



E 
— I 
T 







m 





UE,T) + eV[-^ 



—AT-eV 
T 



Putting this expression into the equations for the electrical and thermal currents, Ean. H1.4fl and Ean. (|1.8|) . one obtains 
the two transport equations 
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These equations are equivalent to the usual form for the transport equations for a metal: 

I = GV + 1-1 AT 
^QV + kAT 

with the linear-response thermoelectric coefficients defined by 
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(1.10a) 
(1.10b) 



(1.11a) 
(1.11b) 



(1.12a) 
(1.12b) 
(1.12c) 
(1.12d) 



If we approximate the derivative of the Fermi function by a 6 function at low temperatures, as we did for the electrical 
conductance G, we see that all other coefficients vanish. In order to obtain a finite value, we use the Sommerfeld 
expansion for the integrals^^ 



= $(0) + 
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dE^ 



E=0 



Using the fact that N'{Q) — N{0)/2Ep, we obtain the following expressions for the last three coefficients 

V 

c 



Ti'^klT eAD N{0) 

6 L 
T^'^ikeTf eAD N{0) 



6 L Ef 

TT^klTAD , , 
^^^-A^(O). 



(1.13) 

(1.14a) 
(1.14b) 
(1.14c) 



Experimentally, the quantities often measured are the thermopower S and the thermal conductance G^. The ther- 
mopower is defined as the thermal voltage generated by a temperature differential AT, under the condition that no 
electrical current flows through the wire (7 = 0). Putting this condition into the first transport equation, Ean. (|1.10a|l . 
we obtain 



V 7/ _ TT^ ks ksT 
AT ^ ^G ^ ~Q~~E^ 



(1.15) 



Note that since the expression for S contains a factor ksT/Ep, the thermopower of a typical normal metal is very 
small. 



4 



The thermal conductance is defined as the ratio of the thermal current to the temperature differential AT, 
under the same condition of no electrical current flow (J = 0). From the equations above 

G^^^^K + CS (1.16) 

For typical metals, the second term is much smaller than the first, and is usually ignored, so that G^ w k. If we take 
the ratio of the electrical to the thermal conductance, we obtain 

Consequently, one finds that the Wiedemann-Franz Law holds, even though the scattering lengths for momentum 
relaxation {£) and energy relaxation (Lin) are quite different. This is because the reservoirs are assumed to be perfect, 
so that any electron entering a reservoir immediately equilibrates with the other electrons in the reservoir. 

Before we go on to discuss normal metals in contact with superconductors, it is worthwhile to review some important 
assumptions in the calculations above. First, in our calculations, we have assumed that the diffusion coefficient D 
is a constant independent of the energy E and position x. When coherence effects are important, as in the case of 
the proximity effect, the diffusion coefficient becomes a function of both these parameters, D{E, x). The differential 
equation for the distribution function, Ean. (|l.ll) is then modified. D(E,x) itself in general is determined by the 
distribution function f{E,x), so one must solve a set of coupled differential equations to obtain a solution. This is 
difficult to do analytically in all but the simplest of cases, and more often must be done numerically. The remainder 
of this chapter will be devoted in great part to deriving the appropriate differential equations for the distribution 
functions and diffusion coefficient in the case of normal metals in contact with superconductors, using the quasiclassical 
equations for superconductivity. 

Second, apart from the electrical conductance G, the thermal coefficients derived above would all vanish if the 
density of states at the Fermi energy N{E), were assumed to be constant, so that it could be taken out of the 
integrals. The small variation in the density of states at the Fermi energy is responsible for finite (but small) values 
of the off-diagonal transport coefficients 77 and C. For example, the thermopower S, which is non-zero only if there 
is an asymmetry between the properties of particles and holes near the Fermi energy, vanishes if N{E) is assumed 
constant. The small difference in N{E) above and below Ep gives rise to the small but finite thermopowers of typical 
normal metals. The conventional quasiclassical theory of superconductivity assumes particle-hole symmetry a priori 
in the definition of the quasiclassical Green's functions, in that N{E) is assumed constant at Ep. Consequently, 
thermoelectric effects cannot be calculated in the conventional quasiclassical approximation; an extension of the 
theory is required. 

Finally, we have been discussing here currents and conductances, rather than current densities and conductivities. 
These are the more relevant quantitites, since we will be discussing mesoscopic samples in which the measured 
properties are properties of the sample as a whole. This will be particularly important for the proximity effect, where 
long range phase coherence means that non-local effects are important. 

2. THE KELDYSH GREEN'S FUNCTIONS 

The starting point for developing the quantum analog to the classical Boltzmann transport equation is the Keldysh 
diagrammatic technique. We shall begin our discussion of the Keldysh technique in the notation of Lifshitz and 
Pitaevskiij2^ As in the equilibrium case, we define a non-equilibrium Green's function 



G,,,,(Xi,X2) = -I (n |tv3,,(Xi)V-+ (X2)| n) (2.1) 



Here Xi and X2 take into account the three spatial coordinates (denoted by fi and f2 respectively), and the time 
coordinate (ti and ^2)- The difference between the non-equilibrium Keldysh Green's function above and the usual 
equilibrium Green's function is that the average is taken over any quantum state |n > of the system, rather than just 
the ground state |0 >. The time ordering operator T has the effect 
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'I (n 
+i (n 



Vi.,(Xi)^+(X2) 

^+(X2)vi.,(Xi) 



n) ii ti > t2, 
n) it t2 > ti 



for fermion operators (the case of interest here) . 
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For simplicity, we label the spin indices and the coordinate indices by numbers. We shall define a set of Green's 
functions that will be useful later: 



G^^ ^-i(Ti;,i;+^, (2.3a) 

G?2^ = z(f^+^2), (2.3b) 

G?2'' = -i(V'iV'2+), (2.3c) 

G^^2"=z(4+^i\. (2.3d) 



The first Green's function is just the one defined above, in this new notation. The second Green's function is similar 
in definition to G"^, except that it is defined with operator T instead of T, which orders the operators in reverse 
chronological order. The last two Green's functions are defined without any time-ordering operators. The four Green's 
functions so defined are not linearly independent, but are related by linear equations of the form 

G"" + G''^ = G"^ + Gf^". (2.4) 

The retarded and advanced Green's functions G^ and G^ can be defined as in the equilibrium case 



-i(4>i-4>t + '4'2'4'i'j ifii>i2, 



and 



G?,^< ^V'"-' ■ ^^^V - (2.5) 

1 iit2>h 



4 1" if ti > t2, , , 




and are related by 



G^2 = (G?i)*. (2.7) 
G^ and G^ can be written in terms of the Kcldysh Green's functions defined earlier as 

qR ^ Qca _ QO^f} ^ QPa (jf3t3 (2.8a) 

(jA ^ _ qPo. ^ Qat3 _ Qf3p (2.8b) 

G"" satisfies the equation of motion 

Go7g(°^""=<5(Xi-X2) (2.9) 
where G^^ is the differential operator (in the free electron approximation)'*^ 



Go-^4 + j!+. (2.10) 



and the second subscript to Gq ""^ denotes that the differentials in Ean. l|2.10|l are with respect to coordinates corre- 
sponding to this subscript. The argument of the delta function Ean. (|2.9l) includes space, time and spin coordinates, 
and the notation G^^"*"" (with a superscript 0) signifies a Green's function for an ideal gas. 

The i5 function in Ean. (|2.9|) arises from the discontinuities in G"" at ti — t2- G^ and G^ have similar discontinuities, 
and obey a similar equation. G^'^^^^ has a discontinuity of the opposite sign, and hence obeys the equation 

GoiGf^^^ = -S{X^ - X2) (2.11) 
G"^ and G^°' have no discontinuities at ti ~ t2, and hence obey the equations 

Goi'g(°/" = (2.12a) 
G^nGft" = (2.12b) 



FIG. 2: Diagrams corresponding to the first order corrections to the Keldysh Green's function G\2 in the presence of an 
external potential U , which is represented by a dashed line. 



a 



a 
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FIG. 3: Dyson's equation in diagrammatic form for the Keldysh Green's function G12 ■ The thick line represents the exact 
Green's function G12*, the thin line represents the Green's function for an ideal gas 0^2°"°", and the ellipse represents the 
self-energy E. 

The diagram technique for Keldysh Green's functions is similar to that for equilibrium Green's functions, except 
that one needs to sum over the internal indices a and (}, with a corresponding increase in the number of diagrams. For 
example, Fig.l shows two diagrams corresponding to the first order corrections to in the presence of an external 
potential, which is represented by a dashed line. 

In the same manner, Dyson's equation for the Green's function can be represented as shown in Fig. 2, where the 
ellipse represents the self-energy function E"*^. In order to make the notation more compact and tractable, it is useful 
to introduce a matrix Green's function and corresponding self-energy matrix: 



(2.13) 

Dyson's equation can then be written in matrix form as 

G12 = G?2 + / G?4S43G32 d^X^ d^Xi (2.14) 



where the usual rules of matrix multiplication are used. 

As we saw earlier, the components of the Green's function matrix are not independent, but are linearly related. 
One can therefore perform a transformation to set one of the matrix components to zero. There are many ways 
to do this; the one we shall use here is the one employed in most recent literature on the quasiclassical theory of 
superconductivity»ii The matrices above are first rotated in Keldysh space using the transformation G —> t^G, where 
the r matrices are identical in form to the Pauli spin matrices 

and then transformed G — > QGQ\ where the matrix Q is given by 

Q=i=(r"-*r2). (2.16) 

The resulting matrices have the form 
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where the retarded and advanced Green's functions and have been defined above, and the new Keldysh Green's 
function G^ is given by 

qK ^ Qaa Q^f) ^ Qaf) ^ QPa ^ (2.18) 

E^, E"^ and are defined in a similar manner. In what follows, we shall use the 'check' to denote the Keldysh 
matrices. In order to avoid writing integrals over the space and time coordinates, we introduce the binary operator 
(g), which has the effect of integrating over the free space and time coordinates and performing matrix multiplication 
when applied between two Keldysh matrices. Thus, the Dyson equation for the Keldysh Green's function written 
above can be represented as 

G^& + & ®t®G (2.19) 

& obeys the equation of motion 

G^^Gl2 = KX1-X2). (2.20) 
Operating G^i on the Dyson equation Ean (|2.19(l from the left, we obtain 

Ggi^G ^ 5{Xi -X2) + t®G. (2.21) 

The conjugate equation is^ 

GGq/ = 5{Xi -X2) + G®t. (2.22) 

Subtracting the two equations, we obtain 

(Gg-i ®G-G® Go 1) - {t®G -G®t) = (2.23) 
where we have suppressed the subscript '1' of Gq^ . This can be written as a commutator 

[Go^eGl-peG] = 

or 

[Go^eG] = [EGG] (2.24) 
where the G operator defines the commutator of two operators [A G -B] = A® B ~ B ® A. 

The mixed or Wigner representation 

The Green's function G12 oscillates rapidly with the difference r2 - n., on the scale of the inverse Fermi wave 
vector kp^. In the physical systems to be discussed, we are interested in variations on much longer length scales. We 
therefore perform a transformation to center-of-mass coordinates, R, T, and difference coordinates, r, t, defined by 
the equations^i 

fi= R- rl2, f2 = R + rl2 (2.25a) 
ti^T-tl2, t2^T + t/2 (2.25b) 

and then define a Fourier transform of G12 with respect to the variables r and t: 

G{R, T; f,t)^ J e~'^'e'P-^ G(i?, T; p, E) dpdE. (2.26) 

In Ean. (|2.24|l . the part of the commutator involving Gp ^ can be written as 

d d \ 1 



G"/ G G - G G G02 



G (2.27) 



^dti dt2 J 2m 

using Ean. H2.10|l in the free electron approximation. Transforming to the coordinates, R^T,f,t^ this can be written 

d 1 



^a^ + -V.--V^~)G (2.28) 
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From Ean. (|2.18() . the Keldysh component of the Green's function can be written as the sum of the Green's 
functions G"'^ and G^". Keeping in mind the definitions of these functions given in Eans (|2.3|l . we can define a 
non-equihbrium distribution function /(i?, T, r), which is related to the function G^" by 



f{R, T,r)^ Nj+{R + r/2, T + t/2)i;{R - r/2, T-t/2) = -iG^'^iR, T, r, t)t=o- (2.29) 
\ / t=o 

Similarly 

U{R + r/2, T + t/2)^+{R - r/2, T - t/2) = 1 - f{R, T,r) ^ iG'^''{R, T, f, t)t=o- (2.30) 
Subtracting the first equation from the second, we obtain 

l-2/(i?,T,r) = /i(i?,r,rO =z(G"''(i?,T,f,t)t=o + G''"(i?,T,r,t)t=o) =«G^(i?,T,f,i)t=o, (2.31) 

where we have defined a new distribution function /i(i?, T, r). In terms of the mixed Fourier transform, Ean. H2.26|l . 
this can be written as 

h{R,T,r) = iG^{R,T,f,t)t^a = 77- / e'P'^ G'^ {R,T;p, E) dpdE. (2.32) 

27r J 

Taking the Keldysh component of Ean. (|2.27|l at i = 0, we obtain 

i'-^^^^'-^^) I dP^^'''G"iR^T;p,E), (2.33) 
or, in terms of the Fourier components with respect to r, corresponding to momentum p, 

dh{R,T,p) , p 



dT 



m 



Vj^h{R,T,p) (2.34) 



where the Fourier transform of the Wigner distribution function h{R,T,p) is given by 



h{R, T,p)^—jG'' (R, T- p, E) dE. (2.35) 

In equilibrium, h is given by 

hoiep) = 1 - 2/o(ep) = tanh(epV2A:Br). (2.36) 

Ean H2.34|) . which is the Keldysh component of the left hand side of Ean H2.24() . has the form of the one side of the 
classical Boltzmann equation for the distribution function. (Using the definition of the function h(R,T,p), this can 
also be written in a more conventional form in terms of f{R, T,p).) The right hand side of the Keldysh component 
of Ean. (|2.24l) must therefore correspond to the collision terms. The right hand side of the Keldysh component can 
be written as 2(S^"G"'^ — E^^'G^"). Taking the limit at < = 0, and writing in terms of the distribution function 
f{R, T,p) using Ean. H2.28|) and Ean. H2.29|) . this Keldysh component of the right hand side of Ean l2.24|) can be written 
as^ 

-2i[j:^'"{R,T,p){l-f{R,T,p)) + j:"f^{R,T,p)f{R,T,p}]. (2.37) 

The first term in Ean. l2.37|l with the factor (1 — f{R, T,p)) has the usual form of a scattering-in term, corresponding 
to the gain of particles, while the second term has the form for a scattering-out term, corresponding to a loss of 
particles. Consequently, we see that the Keldysh component of the right-left subtracted Dyson equation gives the 
transport equation for the distribution function. From the diagonal components components of the same equation, 
one can obtain solutions for the other components of G and S. More typically, the scattering terms on the right hand 
side of the Boltzmann equation make it difficult to solve, and some approximations must be employed. If the variation 
of the system with the center-of-mass coordinates T and R is small, then one can expand the Green's functions and 
self-energies, which are functions of R, f, T, t about R and T in a Taylor's series to first order in r and t. This is the 
gradient expansion discussed by Kadanoff and Baym.^^ and Larkin and Ovchinnikov*, and we shall return to it at the 
end of this section. 
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Instead of taking the difference of Ean H2.21() and its conjugate equation, Ean (|2.22(l . we take the sum, we obtain the 
equation 

[Go 1 ® G] = 2(5(Xi - X2) + [S e G] (2.38) 

where the operator © defines the Keldysh anti-commutator, in the same way as the operator G defines the Keldysh 
commutator. The left hand side of Ean. l2.38|l can be written as 

G-i ® G + G ® G-^ = 

Transforming to the mixed representation, we obtain 

GiR,T,f,t) (2.40) 

Now, the assumption we are making is that the variations of G on the scale of R are much slower than the variations 
on the scale of r. Hence, the the terms in equation above involving derivatives with R contribute much less than the 
those with f, and can be neglected in this approximation. If we consider the equation for Go, for which the terms 
involving S on the right hand side of Ean. l2.38() are 0, then we obtain, after transforming to Fourier components'^ 

{E^ej;)Go{R,T,p,E)^l (2.41) 



d 


- ' ) 


1 


dti 


dt2j 


2m 



G. 



(2.39) 



44 



or 



Go(i?,T,p,i?) = — -, (2.42) 

where 

ep-=- h^vf{p-Pf)- (2.43) 

2m 

So far, we have assumed a free-electron model. If there is a slowly varying potential U{R, T), the equation above can 
be modified to 

Go{R,T,p,E) = ^ (2.44) 

(i?-ep-C/(i?,T)) 

This equation has form of a Green's function for a free particle, but in a slowly varying potential U{R,T). The 
operator Gq ^ in the mixed representation can therefore be written as 

G^^ ^{E-ep--U{R,T)). (2.45) 
More accurately, one transforms Ean. (|2.1()(l using the following equations 

(2.46a) 

(2.46b) 
(2.46c) 



d 


1 d 
~ 2df 


±1 




Vi,2 


1 

-2^« 






V 1,2 


4 


± Vf • - 





Again, assuming the functions in the mixed representation are slowly varying functions of R, we can ignore the second 
derivatives with respect to R, to obtain 
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In most applications of interest here, we also ignore the slow T dependence. Adding the potential U{R,T), and 
writing in terms of the partial Fourier transform with respect to p and E, we obtain 



G^^ ^E-e,j+-VF-Vj^ + fi-U{R,T), 



(2.48) 



where we have replaced p/m by vp, since the important region of interest is at the Fermi surface. Note that we still 
keep account of the direction of the momentum. 

To conclude this section, we derive some expressions for the convolution of two operators in the mixed representation, 
the gradient expansion discussed above. Consider the convolution of two operators defined as an integral over the 
internal space and time coordinates 

A(g)B = dfs dt3v4(fi,f3,ti,t3)B(f3,f2,<3,t2) 



(2.49) 

If we do a partial transform to p, R (but do not transform the time coordinates) , A(E) B can be represented as 

A(^B = J dt3eii^s-^p-^^-^A)A{R,p,ti,t3)B{R,p,t3,t2) (2.50) 

where the superscripts to the derivatives denote that they operate only those functions. If the derivatives are small, 
we need to take only the first order expansion of this expression 



{A(g)B){p,R) 



dU 



1 + _ V - • V- - V- • V - 



Similarly 



{B^A){p,R) = / dt 



1 



R P P R 



A{R,p,ti,t3)B{R,p,t3,t2). 



B{RJM,h)A{R,PMM)- 



(2.51) 



(2.52) 



If we were dealing with functions alone, then the multiplication of the two function A and B is commutative. When 
A and B are matrices, however, they do not in general commute, so we obtain 

[Ae B){p,R) ^ A® B ~ B ® A ^ J dt3[A,B] + ^J dh [{V^A, VpB} - {Vp-A, V^B}] , (2.53) 

where [A, B] notation stands for the commutator of the two functions, and {A, B} stands for the anticommutator. A 
similar equation can be obtained if we transform the times to the mixed representation T, E 



{Aq B){T,E) = A® B - B ® A^ j dr^ [A, B] + ^- j dr^ [{drA, OeB] - {Be A, drB}] , 



When we transform both sets of variables, we obtain 
{Aq B){p,R,T,E) ^ A® B - B® A 



= [A, B] + l- [{{dTA, OeB) - {OeA, OtB}) + ({V^A, VpB) - {VpA, V ^B})] 



(2.54) 



(2.55) 



In most cases, we are interested in stationary situations, where there is no T dependence. In this case, the equation 
above reduces to 

{AeB){p,R,T,E)^A®B-B(DA^[A,B]+^- [{V^A, VpB} - {VpA, V^B}] . (2.56) 
These expressions will be useful in the calculations to follow. 



3. THE QUASICLASSICAL APPROXIMATION 

The non-equilibrium spectral function A is defined in the same way as in the equilibrium caseiii2& 

2tt tt 



(3.1) 
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where 3(G^) denotes the imaginary part of the retarded Green's function. In the equiHbrium case, A defines the 
spectrum of energy levels, for stationary quantum states, it has the form of a sum of ^-functions at each state energy. 
In the quasiparticle approximation, these (5-functions are broadened, but the width T of each state, defining its lifetime, 
is still small compared to its energy. If F is large, the quasiparticle approximation breaks down, and one cannot obtain 
a kinetic equation for a distribution function by integrating over the energy E. However, for most perturbations of 
interest, the self-energies typically have a weak dependence on the magnitude of the momentum, this dependence 
being appreciable only near the Fermi energy. On the other hand, while the Dyson equation has a strong dependence 
on E and Cp-, the subtracted Dyson equation, Ean. H2.24|l . from which we will obtain the equations of motion, has a 
very weak dependence on both E and e^-. In this case, it is possible to average over the particle energy to eliminate the 
dependence on the magnitude of the momentum, but keep the dependence on the direction of the momentum. Hence, 
one can think of replacing the Green's functions and self-energies by their values on the Fermi surface, multiplied by 
a (5-function in the form 

G{R,p,hM) 5{ep)g{R,pMM) (3.2) 
To this end, we define the so-called quasiclassical Green's function 

g{R,p,ti,t2) ^ - dej;G{R,p,ti,t2). (3.3) 



TT 



Care must be taken in performing this integral, since the integrand falls off only as l/ep for large ep. To avoid this, 
one can introduce a cut-off in the integral, as done by Serene and Rainer,^ or following Eilenberger^ use a special 
contour for integration. 

We would like to obtain an equation of motion for the quasiclassical Green's functions. If we obtain an equation of 
motion by operating Gq^ in the form of Ean H2.47() on Dyson's equation, and then integrating over ep, there are terms 
in the equation that will have large contributions. In order to eliminate these large terms, we start from the left-right 
subtracted equation of motion, Ean. (|2.24() . in which the troublesome terms are cancelled, and then integrate over e^-. 
The terms on the right hand side of Eqn. (|2.24|l are of the form 

S(g)G = J dt3dxlt{xi,ti,X3,t3)G{x3,t3,X2,t2) (3.4) 

We first Fourier transform this term with respect to p, which takes care of the integral over . We then average the 
resulting equation with respect to e^:. The assumption here is that only G has a strong dependence on the momentum 
p, so the result of this averaging is a term of the form 

J dt3t{R,p,ti,t3)g{R,p,t3,t2) (3.5) 

Now, to complete the transformation, the self-energy E, which is a functional of the Green's functions G, must become 
a functional only of the quasiclassical Green's functions g, S[G] — > (T[g]. With this final change, Ean. H2.24() for the 
quasiclassical Green's functions becomes 

[{9o'~'^)°9]^0 (3.6) 

The o operator in the commutator [A o B] involves an integral over the internal time coordinates in addition to the 

usual matrix multiplication for Keldysh matrices. If we transform the time coordinates as well, this integral can 
be removed, as shown at the end of the last section fEan l2.55l) . In this case, the commutator becomes a simple 
commutator (but involving matrix multiplication of the Keldysh matrices) . 

It remains to express the physical quantities of interest in terms of the Green's functions. The particle density is 
given by 

p{l) = -2^Gfr (3.7) 
and the current density in the absence of external fields is given by 

j(l)^--(Vi-V2)Gf2"|2=i (3.8) 
m 

From the definitions of G^, G^ and G^ in terms of G"'', G^'^ and G^^^, we can write the function G^" as 
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and depend on the equilibrium properties of the system, and so do not contribute to the current on nonequi- 
Ubrium density. Consequently, these terms can be dropped in the expression for the particle density and current. 
Writing in terms of the quasiclassical Green's functions, we can obtain expressions for the density and the current in 
the mixed representation. Consider the expression for charge density 



dp{R,T) = -ieG^ {R,T,f= 0,t = 0), 



(3.10) 



where we use 6p instead of p to emphasize that this does not include the equilibrium contributions. Expanding G^^ 
in terms of Fourier components in momentum space p 



8p{R,T)^ 



j 4^.^'''G^iR^T,p,t) 



{2n) 

dE 
2^ 



r=0,t=0 



d^p 
(2^ 



G''iR,T,p,E) 



(3.11) 



where we have used the fact that setting t = in G^^ (t) is equivalent to integrating G^ [E) over dE/2TT. 

The equivalent expression for the current density can also be written in terms of the quasiclassical Green's function. 
From Ean H3.8|l and Ean l|3.9|l . we have, in the mixed representation 



j{R,T) = -— (Vi-V2)G«(i?,r,f,p) 
Zm L 



The integral over the momentum p can be rewritten as 

_fp_ 

(27r)3 



No 



dep 



dflp 
An 



(3.12) 



(3.13) 



where iVo is the density of states at the Fermi energy, and dflp is an element of solid angle in momemtum space p. 
Writing Vi - V2 = 2Wp from Ean. H2.46b|l . we obtain 



e 

m 



d^p 
(2^ 



3'P-^G^iR,T,p,t) 



(3.14) 



Operating Vf on the exponential within the integral gives ip. In the quasiclassical approximation, we assume that the 
major contribution comes from near the Fermi surface, so that p/m = vp, the Fermi velocity. With this approximation, 
we obtain 



- , I dE 



d^p 
(2^ 



vfG^{R, T,p,E). 



If we use the definition of the distribution function h given by Ean. (|2.31|l . we obtain 

d^p 



jiR,T) = -e 



VFh{R,T,p) 



(3.15) 



(3.16) 



which is the expected classical form for the current. Recalling the definition of the quasiclassical Green's function, 
Ean H3.3|l . we obtain 



-eNf 



j dE j ^vpgK(E,p,R,T) 



(3.17) 



where we have used Ean (|3.3|l and Ean. (|3.13|l . Note that writing these expressions in terms of the quasiclassical 
Green's functions necessitates reversing the order of integration over E and e^. 

The transformation, Ean. H3.13|l . assumes that the density of states is constant at the Fermi energy, and hence 
assumes that particle-hole symmetry holds. Hence, within this approximation, we will not be able to obtain any 
results on physical phenomenon that depend on particle-hole asymmetry, in particular, thermoelectric effects. This 
should be contrasted with the derivation of the current in the diffusive limit that we performed in the introduction, 
where the energy dependence of the density of states was taken into account explicitly. 

So far, we have ignored the effects of external fields and potentials. In particular, the effect of a magnetic field is 
of interest. The magnetic field is introduced in the form of a vector potential A{R,T). Here, we shall consider only 
time independent fields. A{R) is introduced by making the change 



(3.18) 
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in all equations involving the spatial derivative. For example, the operator Gqi is now written 

G-=.|^ + g+M-e0. (3.19) 

In the equation above, we have also added a term e(f), corresponding to the presence of a scalar potential. 

In our formulation, we would like all observable quantities to be invariant under a gauge transformation of the 
vector potential 

A + Vx{r), (3.20) 

which also transforms the potential 



|. (3.21) 



Ean. (|3.14|l for the electrical current is then modified to 

2e 



m 



(2^ 



[Vp-ieA{R)] / e'P"G^"(i?,T,p,i) 



(3.22) 



where we have written the current in terms of G^"' rather than . It follows that Eon. H3. 15(1 can be written as 

J(i?,T) ^ -.2e / f / If, ^^.G^-iR,TJ.E) + ^I(i?) / f / |f, (i?, T, i.) . (3.23) 

The second term in the equation above is called the diamagnetic term. It is cancelled by a contribution from the first 
term arising from energies far from the Fermi energy, which are not taken into account in the quasiclassical Green's 
fimction as defined by Ean. (|3.3|) . — Note that such a cancellation does not occur in the case of a superconductor, and 
the second term gives rise to the supercurrent, which is proportional to the vector potential and the phase gradient 
according to Fan. ((3. 21(1 . With this high energy contribution cancelled by the diamagnetic term, the equation above 
transforms into the expression Fan. ((3. 171) for the electrical current, written now in terms of the Keldysh component 
of the quasiclassical Green's function. 

Such a cancellation does not occur in transforming 5p in terms of the quasiclassical Green's functions, and the 
contribution of the integrals in Fqn. ((3.11() must be explicitly calculated]— The result that was obtained by Eliashberg^ 
can be written as 

Sp{R,T)^^^ J dEj -^g''{E,p,R,T)-2e^Nocj){R,T), (3.24) 

where 4'{R, T) is now the scalar electrochemical potential. This expression can also be obtained by invoking gauge 
invariance arguments. Since we must preserve charge neutrality in the system, 5p = 0, so that (j){R,T) is given by 

(/.(i?,r) ^-^jdE j ^g^{E,p,R,T) (3.25) 

For completeness, we also write expressions for the thermal current and the density of states. The expression for 
the thermal current in terms of the Keldysh Green's function is given by 

lh{R,T)^i J^J ^ EvfG''{R,T,p,E), (3.26) 

and the corresponding form in terms of the quasiclassical Green's functions 

jthiR,T) = iiVo JdEj ^EvFg''iE,p,R,T). (3.27) 

As in the equilibrium case, the density of states is given directly by the spectral function A defined in Eon. 1(3. 1(1 . 
now expressed in the mixed representation 

N{E, R, T) = A{E, R, T) = --'^G"{E, R, T) 

TT 



= / J2^G'\EJ,R,T). (3.28) 
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Written in terms of the quasiclassical Green's function, this becomes 

N{E, R, T) = iVoSR j ^g'^iR: E, T), (3.29) 
where the notations 3? and 3 stand for the real and imaginary components respectively. 

4. NON-EQUILIBRIUM GREEN'S FUNCTIONS FOR SUPERCONDUCTING SYSTEMS 

The formalism that we have developed for non-equilibrium Green's functions for normal systems carries over into the 
superconducting case, except that the Green's functions for the superconducting case are more complicated. In order 
to show how the Green's functions are defined, we start by expressing the Green's functions in the superconducting 
case in terms of field operators in the Nambu-Gorkov formalism. In this formalism, the field operators can be written 
as two-component column matrices 




*i = (4.1) 



where the up and down arrows refer to the spin indices, and the numbers now refer only to the time and space 
coordinates. The Hermitian adjoint of this operator can be written as a two-component row matrix 

= (^+ V^u) • (4.2) 
The multiplication of two such operators is defined as the tensorial product of the two matrices. For example 

^ ( f^fi] (4.3) 

It is natural then to define the Green's functions in terms of these products. For example, the natural definition of 
the Green's function corresponding to G"", Ean. H2.3al) . would be 

However, we would like to keep the same form for the equations of motion and Dyson's equation as for the normal 
case, except, of course, the quantities will be matrices in Nambu-Gorkov space. To see if the definition above fits this 
requirement, let us operate id/dti on the Green's function for an ideal gas, as defined above. With the definition 



_5_^(0)oa_, / (^V^IT^'t) fS{X,-X2) 



we have 

This does not have the required form of the analogous equation for the normal case, Ean. H2.9|l . To remedy this, we 
define our Green's functions with an extra factor of t^. The definitions corresponding to Eans. (|2.3|l are then 



= -zr' = -zr^ I ) ^ ( ) ^ 77 I , (4.7c) 




G?2" = -ZT^ (T^,^+ ) = -zt' I )^ \- ( \- 'r{ \ , (4.7a) 

Gfl = -zr' (f4>,i>t ) = -zr' I )^ i )^ ri i , (4.7b) 



15 



and 



G 



13a 
12 



' V'2i'0]*, 



(4.7d) 



where the 'tilde' over the Green's functions denotes that they are matrices in Nambu-Gorkov space. The operator 
corresponding to Ean. H2.9|l is then defined as 



"1 dti 2m 



(4.8) 



where it is understood that any 'scalar' quantities in the equation above and in what follows are mulitiplied by 
the identity matrix r°. With these modifications, all the equations derived for the Keldysh Green's functions for 
the normal case can be carried over directly to the superconducting case. In particular, the advanced and retarded 
Green's functions and G^ are defined in terms of the Green's functions in Ean. H4.7|) in the same manner as before. 
The Keldysh matrices corresponding to Ean. (|2.17|l are then 



G 



qR qK 

G^ 



E-^ 



(4.9) 



Since each element of these matrices are themselves 2x2 matrices, the resulting Keldysh matrices for superconductors 
are 4x4 'supermatrices,' and we shall denote them by a 'hat' symbol ( " ). The equation of motion equivalent to 
Ean. (|T^ is 



[Go'eG] = [seG]. 



(4.10) 



Before we go any further, we need to specify the self-energy S. At the temperatures of interest in the experiments 
on proximity systems, the important self-energy terms are due to electron-phonon scattering (Se_p), and electron 
impurity scattering (Eimp). For conventional superconductors, the elastic component (S^_p -f S^_p) of the electron- 
phonon contribution is the one that leads to the coupling between superconducting electrons^ (the other contributions 
of electron-phonon scattering will be ignored, under the assumption that we are at low enough temperatures so that 
inelastic electron-phonon scattering can be ignored). Perhaps a simpler way of dealing with this component of 
the electron-phonon interaction is to start directly with the Gorkov equations of motion^^ for the Green's function 
Ean. (|4.7a|l . which we write in the form (ignoring impurity scattering for the moment) 



'9*1 



2in 



A 



2m 



5{Xi-X2), 



with the pair potential A defined as^ 



A = A liniit < Til)i^i}j2i >= i\ lim^it [G?2"]i2 



(4.11) 



(4.12) 



in terms of the upper left component of the Green's function G"^, which is frequently called the anomalous Green's 
function (or pair amplitude), and denoted by F . Here A is the coupling constant. For a uniform bulk superconductor 
A is real. For a normal metal, A vanishes, and hence, although the pair amplitude in a normal metal may be finite, the 
pair potential vanishes. Ean. ()4.12fl defines a self-consistent equation for the pair potential A; it is defined in terms of 
the anomalous Green's function i^, which in turn is determined by an equation that depends on A. Note that unlike 



?r2'']i2 or [Gf^"] 

Still ignoring impurity scattering, the equation of motion Ean. H4.1U|l can then be written in compact form as 



[G"^]ii, [G"^]i2 is continuous at ^2 = ^i, so that A can also be written in terms of [G^j ]i2 or 

[G?2"]l2 ati2=ii. 



IT 



d_ 



2rn 



M + A eG 



0. 



(4.13) 



Here A represents a 4x4 matrix 



A 
A 



(4.14) 
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where 

with A defined by Ean. H4.12ll . and 

- (^n ] ■ (4.16) 



Before we move on to making the quasiclassical approximation, it is useful to obtain expressions for measurable 
quantities in terms of the Green's functions defined in this section. Following Eqn. (|2.32(l , one can define a distribution 
function by averaging the Keldysh component of the Green's function defined in Ean. H4.9|l 



hiR,T,p)^—JdEG''iR,T;p,E). (4.17) 

However, since is a 2x2 matrix, h is also a 2x2 matrix, so that its interpretation as a simple distribution function (in 
the flavor of Ean. (|2.32|l for the equilibrium case, for example) is not immediately clear. A more physical interpretation 
can be obtained by diagonalizing the Gorkov equations given in Ean. H4.11|l i^ The matrix on the left hand side of 
this equation can be diagonalized by a unitary transformation; this transformation, of course, is just the Bogoliubov- 
Valatin transformation, which results in a diagonal 'energy' matrix with eigenvalues 

El = e}+A\ (4.18) 
The same transformation also diagonalizes the equilibrium distribution matrix ho, which now has the form 

, ^/1-2/ot(£;,^ \ 

2foi{Ep)-l)- ^^■^■^> 

The top-left component is for electron-like excitations, and the bottom-right component for hole-like excitations. This 
suggests that for our combined Nambu-Gorkov-Keldysh Green's functions, the equations for the electric current and 
thermal current should be modiflcd to 



le 



dE f d^p 



j{R,T) = -- J ^vfTt t'G''{R,T,p,E) , (4.20) 



and 



dE f d^p 



jMT) ^--j^j ^EvpTr [g''{R,T,p,E)\ . (4.21) 

Taking the trace of the matrix takes the contributions of both electrons and holes. Note that the expression for 
electrical current includes a factor of in the argument of the trace, while the thermal current does not. 



5. QUASICLASSICAL SUPERCONDUCTING GREEN'S FUNCTIONS 



In principle, the properties of the superconducting system may be calculated starting from Gorkov's equations. 
In practice, however, such calculations are difficult in all but the simplest cases. The problem is that Gorkov's 
equations contain information at length scales much finer than those of interest. The way around this is to make the 
quasiclassical approximation as we did for the case of the normal Green's functions, which we proceed to do below. 

Taking into account impurity scattering, the equation of motion for the superconducting Green's function can now 
be written as 



2m 



A-E. 



tmp 



eG 



= 0, 



(5.1) 



where the impurity self-energy S^mp includes contributions from both spin-flip and spin-independent elastic scattering. 
The elastic contribution to the self-energy can be written as 



d^p' 

(27r)3 



vip-p') G{p') 



(5.2) 
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Here V is the impurity potential, and Ni the number of impurities per unit volume. We assume that v{p) is independent 
of the magnitude of p, so that 

Eo(p) = N^No I dep ^ Hp-p')]"" G{ff). (5.3) 
Defining the elastic scattering rate l/r in the Born approximation by 

- = 2nN,No^\v{p-p')\\ (5.4) 



we can write 

° 2ttt 



depG{p) (5.5) 



Here Giji) is a function of the magnitude of p alone. Similarly, for the contribution from spin-flip scattering, one 
obtains 

tsf = -^ I dtpf^G{p)T\ (5.6) 

sf J 

From EQn. H4.9(l . one has 



if^^ + -^+fi + A-±o-±sf]eG 
oti 2m ' 



0. (5.7) 



Converting the left hand side of this equation to relative coordinates using Ean. 1)2.46(1 (neglecting any derivatives with 
respect to the center-of-mass time T), and using only the first term of Ean. H2.56l) to lowest order, we obtain 

[t^E + a, G] + ivp ■ dgG - [Eo + S^/, G] = 0. (5.8) 

Note that there are no integrations over space or time coordinates in the equation above, but only matrix multiplica- 
tions. Written in terms of the quasiclassical Green's functions, this becomes 

[t^E + A,g\+ ivF ■ d0 - [ao + cr,f,g] = 0, (5.9) 

where 



and 



CTo = (5-10) 
^s} = ~i-^T^,T\ (5.11) 

2Tsf 

where the subscript to the quasiclassical Green's functions denotes that they are averaged over all angles of p. Eqn. (|5.9() 
is Eilenberger's equationji and is the starting point for many calculations on superconducting systems. 

To complete the transformation to quasiclassical Green's functions, we express the equations for physically ob- 
servable quantities in terms of the quasiclassical Green's functions. From Ean. (|4.12|l . the gap parameter A can be 
expressed as 

A^zA[G?f+]i2 = ^A[GfiJi2 (5.12) 



using the definition Ean. l2.18|l of the Keldysh Green's function. In the mixed representation, this can be written as 

= N„-JdE I |i[g''(fl.i,p,E)],, 

JdEj ^^Tr{{r'-,T')g'<lfi.t.p,E}] (5.13) 



18 



From Eans. H4.20() and (|4.21(l . the electrical current and thermal currents are given by 



j{R,T)^ / dE / -r^vpTv T'~g^{R,T,p,E) 



and 



No 



An 



jth{R,T) dE -^EvpTv ~g''{R,T,p,E) 



Aw 



(5.14) 



(5.15) 



One can also define an equation that gives the amount of charge associated with quasiparticle excitations. In a 
superconductor in equilibrium, the number of particles and holes are equal, so the quasiparticle charge should vanish. 
If there is an imbalance in the population of electrons and holes, one can obtain a charge imbalance^'^ usually denoted 
by Q*, which is given by 



Q* 



le 
2" 

eiVo 
4 



dE 
2^ 



dE 



d^p 
(2^ 

da 

Att 



Tr 



^Tr 



G^{R,T,p,E) 
5^(i?,T,p,i?) 



(5.16) 



Q* is essentially the first term on the right hand side of Ean. (|3.24|l . Invoking charge neutrality, this then would result 
in a electrochemical potential given Ean. H3.25() 



cl,{R,T) = -^Q* 



(5.17) 



The normalization condition and the distribution function 

In obtaining the equation of motion of the Green's functions by subtracting the Dyson equation from its conjugate, 
some information was lost regarding the norm of the quasiclassical Green's function g. The norm of g can be obtained 
by the normalization condition for the quasiclassical Green's function obtained by Eilenberger^ and Larkin and 
Ovchinniko\ii 

99 (5.18) 

where (4x4) matrix multiplication is implied. This normalization condition can be obtained by explicit calculation 
for a bulk system in equilibrium. Furthermore, it can be shown that the normalization condition is consistent with 
Eilenberger's equation Ean. H5.9() . 

Ean. H5.18|) is equivalent to the three (2x2) matrix equations 

= (5.19a) 

(5.19b) 

and 

5^5^ = 0. (5.19c) 

As in the case of the i?-integrated Green's functions, we would like to obtain an equation of motion for a distribution 
function, eqiuvalent to the quantum Boltzmann equation derived earlier. We introduce such a distribution function 
h by the ansataii 

g^=g^h-hg^. (5.20) 

This form of the quasiclassical Keldysh Green's function satisfies the normalization condition Ean. H5.19c|l . as can be 
verified by direct substitution 

f ig^'h - hg^) + {f h - h~g^)g^ = (5.21) 
using Eans. H5.19a|l and H5.19b|l . Note that the function h is not uniquely defined; if h is replaced hy h + h', where 

h' = g'^i + S:g^ (5.22) 
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and X is an arbitrary matrix function, the right hand side of Ean. (|5.20() is unchanged.^ This arbitrariness allows us 
some flexibility in choosing the distribution function h. At low frequencies, for example, following Schmid and Schon^ 
and Larkin and Ovchinnikov)^ we may choose h to be diagonal in particle-hole space 

h = hLT" + hTT^. (5.23) 

The subscripts refer to the longitudinal (/iz,) and transverse (/it), a terminology introduced by Schmid and Schon 
to refer to changes that are associated with the magnitude (/i^) or phase (hr) of the complex order parameter. In 
equilibrium, /il(£^) = 1 - 2foiE) and hriE) = 0. 

Another possible choice is one introduced by Shelanko'ViSS 

h^hiT°+g^h2. (5.24) 

This representation has the advantage that the equation for the distribution function reduces to a form very similar to 
the Boltzmann equation when the quasiparticle approximation is valid. In what follows, we shall use the representation 
of Schmid and Schon. 

The representation of the Keldysh Green's function in terms of a distribution function allows us to obtain an equation 
of motion for the distribution function from the equation of motion for the Green's function. From Eilenberger's 
equation H5.9|l . we obtain three equations of motion for the three components of the quasiclassical Green's function 

[Et^ + A, 5«] + ivF ■ dj^g^ ~ [d^, g«] = (5.25a) 
[Et^ + A, ~g^] + iVF ■ d^g^ ~ [B^ ,{1^] = (5.25b) 
[Et^ + A, 5^'] + ivF ■ dj^g'' - d'^cj^ - ^^.g^ + + g^a^ = (5.25c) 

where a — + dsf- Substituting g^ from Ean. H5.2()|l into Ean. l|5.25c|l . we obtain 

[Et^ + A, g^h - hg^] + ivp ■ dg{g^h - hg^) - d^{g^h - hg^) - a^g^ + g^a'^ + {g^h - hg^)d^ = 0. (5.26) 

Subtracting from the equation above Ean. H5.25a|) multiplied by h on the right, and adding Ean. lj5.25b|) multiplied by 
h on the left, we obtain 

-g'^Blh] - B[h]~g^ = 0, (5.27) 

where 

B[h] = [Et^ + a, /i] + ct^ - {d^h - hd^) + ivF ■ dj^h. (5.28) 

Eauation (|5.27|l is the required equation of motion for the distribution function h. 

As an example, we shall calculate the Green's functions for the equilibrium case for a bulk superconductor, in the 
limit where a — Q. Eauation H5.25a|l in this limit has the form 

[Et^ +K,g^]=Q (5.29) 

First, if we represent the retarded Green's function by 

- (:■; :S) • ("») 

the normalization condition Ean. H5.19a|) and Ean. (|5.29|) together imply that 

511 = -.922, (5.31a) 

5ii +521312 = 1, (5.31b) 
A* 

521 ==-312^ (5.31c) 



and 



E 

511=512^- (5.31d) 
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Solving these equations, we obtain 

( E A 
A* g 

In taking the square roots in the equation above, a factor of i wiU appear if i?^ < |Ap. This is the expected solution 
from Gorkov's equations. Note that gn is just the normalized BCS density of states for E'^ > |Ap, so that gn 
represents a generalized density of states. 

From the equations above, it is clear that 1721 = ~9i2- To obtain g^, we can use the relation Ean. l)2.7|l . which in 
terms of Nambu-Gorkov matrices reads 

G^^-T^G^'+T^ or 5^ = -rV+T3 (5.33) 

where the ^ symbol denotes the Hermitian conjugate. 




6. THE DIRTY LIMIT: THE USADEL EQUATION 



In systems where elastic impurity scattering is strong, the motion of quasiparticles is not ballistic, but diffusive. 
In this case, the effect of the strong impurity scattering is to randomize the momentum of the quasiparticles, so that 
it makes sense to average the properties of system over the directions of the momentum, keeping in mind that the 
magnitude of the momentum is still its value at the Fermi surface, pp- 

With this in mind, let us consider a system with strong impurity scattering, which we define as a system in which 
the impurity scattering rate l/r, as defined by Ean. H5.4|l . is much larger than any energy in the problem (i?. A). 
In this case, we expand the quasiclassical Green's function to first order in momentum (essentially an expansion in 
spherical harmonics) 

9 = 9s+pgp (6.1) 

where p denotes a vector of unit magnitude in the direction of p, and gs and gp are independent of the direction of p, 
and the subscripts stand for s-wave and p-wave components of g. The assumption is that gp gs- The self-energy is 
expanded in a similar fashion 

(7 = (To + (7s/ ^ (Ts+ P(7p (6.2) 

We would like to calculate the components of a in terms of the components of g, for which we shall need the definition 
of a in terms of g. Ignoring spin-flip scattering for the moment, this relation is given by Ean. (|5.3|l . which can be 
rewritten in terms of the quasiclassical Green's functions in the form 

aoip) = -tTTN.No I ^\v{p-ffg{jf) = -^TriV.iVo J ^\v{p ■ p')f{gs + p' gp,), (6.3) 

using the expansion Ean. (|6.1|) for the Green's function above. In order to do this, we take the dot product of p with 
both sides of the equation above 

p ■ a = p&s + (p ■ p)^p = P^s + d-p. (6.4) 
We perform a similar operation on Ean. (|6.3|l 

-^\vip-p')\''ipgs + {p ■ p')gp'). (6.5) 

Since both p and gs are independent of the integration over dflp' , the first term under the integral in the equation 
above can be written as 

-iirpgsN^No J ^\v{p ■ p')\^ - ~^pgs, (6.6) 

using the definition Ean. H5.4|l of l/r. If we consider non-spin-flip scattering alone, then a can be obtained by equating 
like terms in Eans. H6.4|l and (|6.5|l . If we include spin-flip scattering, one can write 
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as expected from Eans. (|5.10|l and H5.11|l . For the p-wave component, we consider only the contribution from non 
spin-flip impurity scattering, under the assumption that it is much stronger than the spin-flip scattering. The second 
term under the integral in Ean. (|6.5(l can be written as 

-ZTTiV.TVo J ^Hp-p')\^{p-p')gp,^~t7TN,No9pJ ^\v{p ■ p'^l - {1 - p ■ p% (6.8) 

where gp can be taken out of the integral since it is independent of the direction of p. The first term in the square 
bracket can be seen to be related to the elastic scattering rate l/r. The remaining terms can be written in terms of 
the transport time, defined by 

— = 27riV,7Vo / ^\vip.p')\^il-p.p'), (6.9) 

Ttr J 47r 

that is well known in the transport theory of metals. With this definition, cTp can be written as 

^P = -k--—)9p- (6.10) 

2 T Ttr 

By putting Ean. H6.1|l into the normalization condition for the quasiclassical Green's function, and neglecting terms 
quadratic in gp, we also obtain the two equations 

gs9s = 1 (6.11a) 

and 

9s9p + 9p9s = 0- (6.11b) 

We now proceed to expand the Eilenberger equation, Ean. (|5.9f) . in terms of the s and p-wave expansions of g and (t. 
Replacing by vpp, we obtain 

[t^E - A,gs] +p[t^E ~ A,gp] +ivFp- dj^{gs + pgp) - [a-s,gs] ~p[a-s,gp] + [o-p,gs] + p-p[d-p,gp] = 0. (6-12) 

The last term is second order in the small quantity gp, and can be neglected. Collecting the terms that are even in 
p, we obtain 

[t'^E + A, ,9,] + ivpip ■ p)dggp = 0. (6.13) 
Averaging this equation over all directions of p gives 

[T^E + A,gs]+i^d0p^O. (6.14) 

Ignoring spin-dependent scattering for the moment, the terms that are odd in p can be written 

[T'^E + A,gp]+ivFd0,--^[gp,g,]=O, (6.15) 

where we have used EQns. (j6.7|) and H6.10() to write as and ap in terms of gs and gp. If elastic scattering is strong, the 
first term in the equation above can be neglected compared to the third, so we obtain 

ivpd^s H gsQp = 0, (6.16) 

where we have used Ean. lj6.11b(l to simplify the second term. Multiplying this equation by cjs on the left, and using 
Ean. (|6.11a|l . we obtain 

VFgsd0s = — —gp, (6.17) 

or writing gp in terms of g^ 

gp = vpTtrgsd^gs = -Igsd^s (6.18) 
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where we have introduced the elastic scattering length ^ = vpTtr- Putting this into Ean. (|6.15|l . we obtain 

[T^E + A,g,]^i^d0sd0s = O, (6.19) 
or writing this in terms of the diffusion coefficient D (1/3)vf£, and reintroducing the spin-flip scattering term, 

[f'^E + A - asf,9s] - iDd0sd0s = 0. (6.20) 

This is the equation first derived by Usadel^^S and forms the starting point for most discussions on dirty supercon- 
ducting systems. 

In the remainder of our development, we shall neglect the spin-flip scattering term. Writing gs as a matrix 

9s = ?a] , (6.21) 



.05.. 

we can write the matrix equation, Ean. (|6.2Q|l . as three separate equations 

[t^E + A, ] = zDd^{g^d0^), (6.22a) 



[t'E + A, ] = zDd^{~gtdj^~gf), (6.22b) 

and 

[t'E + A, gf] = zDd^ [(5f 5^5f ) + (g^d^f)] . (6.22c) 

As before, the first two equations come from the diagonal components of the Usadel equation, and describe the 
equilibrium properties of the system, while the third equation comes from the off-diagonal or Keldysh component, 
and represents the kinetic equation for the distribution function. These equations are supplemented by the equations 
for measurable quantities corresponding to Eans. H5.14|l . I|5.15(l and (|5.14(l . Expanding g^ in Ean. (|5.14l using Eans. H6.1|l 
and Eons. H5. 14(1 . we have 

J(i?,T) = -'-fJdEj ^vppTr P {g^+p~g^)] = ^ / dETr [r^~gsd^g.f] , (6.23) 

where the angular average over the p-wave component give a factor of (1/3), and we have replaced (1/3)wf^ by the 
diffusion coefficient D. Now 

(g,9^g,)^ = (fff S^fff ) + (g^d^f) (6.24) 
so that Ean. (|6.24|l above can be rewritten as 

j{R, T)^'^j dETr P (5f a^5f + 9^0^^)] . (6.25) 

The equation for jth is obtained in a similar way 

jM T) = ^ j dEETT [gfd^gf + ~gfd^gf] , (6.26) 

The kinetic equation, Ean. H6.22c|l . can be recast in terms of a differential equation for a distribution function /i, in 
the hopes of separating the equilibrium properties of the system, represented by g^ and gf , from the non-equilibrium 
properties of the system, represented by h. To this end, as before, we substitute Ean. H5.20|) for g^ in Ean. (|6.22c|l . 
which then becomes 

[t^E + hrgfh] ~ [t^E + A, h~gf] = iDd^ [{gfid^fYh - h~gfid0f)) + d^h - [gf'd^h~gf)] . (6.27) 
Taking the trace of both sides of the equation above gives 

= iDdjiTT[dj^h + ~gfid^~g^)h - h~gj{d^gf) - ~g^{d^h)gj], (6.28) 
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using Eans. (|6.22a|l and (j6.22b|) . In the linear regime, we can use the diagonal representation of h given by Ean. (|5.23|) 
to obtain 



= 0, 



where we have used the fact that 

Tr{9^(gf gf )} = Tr{(5^(5f )5f + fff )} = 2Tr{gf )} = 0, 

Taking the trace after mulitiplying both sides of Ean. H6.27|l by gives 

Dd^ [d^hrTril - fff r^g.^r^} + /j^Trlr^ {~gf {d^~gf ) - gfiO^f)) } - d^hLTv{g^~g^T'}] 

= I kTr{r3[.9f - ~gf, A]} - 2/.TTr{ A(gf + .9,^)} 



(6.29) 



(6.30) 



(6.31) 



Equations (|6.29|) and (|6.31l) form a set of coupled differential equations for the distribution functions h l and Ht ■ Let 
us define the quantities 
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Tr{r3 {9f{.d^gf)-~g^{d^gf))] 



and 



Then Eans. H6.29|l and H6.31(l can be written in the form 
dii [Mo^id^hL) + Qhr + Maoa^/ir)] = 0, 

[M33(9^/iT) + QhL + Moaa^/iL)] = ^ [/iLTr{r3[gf - ~gf , A]} - 2;iTTr{A(gf + ~gf)] 



(6.32) 

(6.33) 

(6.34a) 
(6.34b) 



These equations are in the form of diffusion equations for the distribution function, more general forms of the diffusion 
equation discussed for the normal case in the introduction. As we shall see, the quantity Q is related to the spectral 
supercurrent in the system, DMij are now the energy and position dependent diffusion coefficients, and Eans. Hfi.34|l 
are essentially continuity equations for the spectral thermal and electric current. In the normal limit, g^ = r^, 
5^ = -T^, and A = 0, so that Moo = M33 = 1, and Q = M03 = M30 = 0. Equations (|6.34|) then reduce to 
Ean. (|l.ll) . as expected. Noting that the term in square brackets in Ean. (|6.22cl) is the same as the term in parenthesis 



(6.35) 



in Ean. H6.25|l and the term in square brackets in Ean. (|6.26|l . the electric current can be written as 

j(i?,r) = eNoD [ dE {A-h^idghr) + QHl + Mo^dghL). 



The first term corresponds to quasiparticle (or dissipative) current, and the second term to the supercurrent. The 
third term, which is proportional to the derivative oi Hl, is associated with an imbalance between particles and holes. 
The thermal current can be written in a similar way 

jthiR, T) = NoD J dE E[Moo{d^hL) + Qhr + M^^od^hr]. (6.36) 

For the charge- imbalance Q* , we note that only the s-wave part of the Keldysh Green's function in the square brackets 
in Ean. l5.16|l survives after angular averaging. Writing gf in the form given by Ean. H5.20|) . with h given by Ean. H5.23|l . 
we have 

Q*=-^JdE hLTT{gf ~ gt} ~ ^ J dE hrTvir^ig^ - g^)} = -eNo J dE hrNiE), (6.37) 
since gf and gf are traceless. Here we have defined the normalized superconducting density of states by 

N{E) = lTT{r^g^-g^)}, (6.38) 
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which reduces to the conventional BCS density of states in the equihbrium case. 

We can also recast the kinetic equations Ean. (|6.34|l in a slightly different form sometimes used in the literature. To 
do this, we subtract Ean. l)6.22b(l from Ean. H6.22a|l . multiply by r^, and take the trace. The result is 



Using the definition of Q, we have 



I 

AD 



TrP[5f-5^A]}. 



(6.39) 



(6.40) 



If we choose a gauge in which A is real, the right hand side of the equation vanishes, so that dj^Q — 0. Clearly, 
d^Q — also for a normal metal (even a proximity-coupled normal metal), where A vanishes. In either case, the 
spectral supercurrent Q is conserved. The right hand side of the equation above multiplied by is the same as the 
second term of Ean. lj6.34b|l . Subtracting Ean. H6.39() multiplied by from Ean. (j6.34b|l . we obtain 



[M33(a^/iT) + ^WghL)] + Qid^hL) = [2hTTT{A{~g^ + ~g^)} 

For real A, when dj^Q = 0, Ean. (|6.34a|l can be written in a similar manner 

dji [Mooidj^hL) + MsoidjihT)] + Qidj^hr) = 
Finally, we can also write the kinetic equations in a form similar to Ean. (|5.27|l 

~g^B[h] - B[h]g^ = 0, 



(6.41) 



(6.42) 



(6.43) 



by performing similar manipulations on Eans. (|6.22|l as were performed on Eans. (|5.25|l . For the diffusive case, B[h] is 
a functional of h now given by 



B[h] = [T^E + A,h] -iD 
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(6.44) 



Boundary conditions for the quasiclassical equations of motion 

Eans. l|6.22|) or their derivatives form a set of coupled differential equations for the quasiclassical Green's functions and 
the distribution function. In order to obtain a solution, however, we need to specify the boundary conditions for the 
Green's functions and the distribution function. To this end, we define the concept of a reservoir, where the Green's 
functions and distribution function have well-defined values. 

For a normal reservoir, the retarded and advanced quasiclassical Green's functions are given by 

9m = 9m = (6-45) 

and for the superconducting case, by Ean. H5.32|) . which we reproduce here 

^« = f I . (6.46) 

with given by Ean. (|5.33|l . 

The equilibrium distribution function h is given by Eqn. (|4.19|l in both normal and superconducting reservoirs (since 
we are dealing only with excitations), where, as we noted earlier the (1, 1) component of the matrix applies to particle- 
like excitations, and the (2, 2) component of the matrix applies to hole-like excitations. From this point of view, the 
Fermi functions in Ean. (|4.19|l are given in terms of the usual equilibrium Fermi function Eqn. 511 by fa^{E) — fo{E) 
and foi{E) = fo^—E).^^ Looking ahead to where we might have a finite voltage V on a, reservoir, the equilibrium 
form of h can then be written as 

/ tanh f 



tanh 



E-eV 
2kBT 
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If we write ho in the form of Ean. l|5.23() the equiHbrium values of and Ht can then be expressed as 

If a finite voltage is put on the superconductor, we will obtain a time evolution of the phase in accordance with the 
Josephson relations. Since we have restricted ourselves here to the static case, we must assume that the voltage on 
the superconducting reservoir = 0. In this case, hr = for the superconducting reservoir. 

For a system with perfect interfaces between the superconducting and normal parts, the boundary conditions defined 
by the equations above are sufficient to solve the differential equations, the implicit assumption being that the Green's 
functions and the distribution functions are continuous across any interface. When the transparency of the interface 
is less than unity, this is no longer true. Zaitse\*0^ derived the boundary conditions for the Green's functions at an 
interface of arbitrary transparency. Kupriyanov and Lukichev^^ simplified these equations for the diffusive case in 
the limit of small barrier transparency. Consider then an interface at x = between two metals, say one a normal 
metal in the half-plane x <0 (labeled by the index 'I') and one a superconductor in the half-plane x > 0, (labeled by 
the index '2'), although it also could be an interface between two different normal metals or superconductors. The 
boundary conditions of Kupriyanov and Lukichev are then 

VFiDigsiid^gsi) = VF2D2gs2{dxgs2), (6.49a) 

9sA9si ^ ^[9si,9s2]- (6.49b) 
2r 

Here dx denotes a derivative in the positive x direction, and r — Rh/RN is a parameter that is nominally the ratio 
of the barrier resistance Rb to the normal wire resistance per unit length Rn/L, but is inversely proportional to the 
transmission t of the interface. The first of the equations is clearly related to the conservation of current across the 
interface. The right hand side of the second equation has been shown to be the first term in an expansion of terms of 
the transmission coefficient and hence, it is only valid for low t. 

The diagonal part of Eqn. (|6.49b|) gives the boundary condition for the Green's functions 

5fi5.5fi = ^ m~9^2 - ~9^29^i] (6.50) 
The off-diagonal part of Ean. lj6.49bp gives the boundary condition for the distribution function 

9?id.9^i + g^id.gfi = ^ [5^15^2 + 9^i~9i " 9^29^1 " 3.^5.1] (6.51) 

If we put in Ean. (|5.2U|l for , with Ean. (|5.23|) for h, and then take as before the trace of the resulting equation, and 
the trace of the equation multiplied by , we will obtain boundary conditions for and hj^ . Noting that the left 
hand side of Ean. (|6.51|l is simply the term in square brackets on the right hand side of Ean. H6.22c|l . we obtain the 
two equations 



and 



2r [MooidghL) + Qhr + ^h^dghT)] = a^, (6.52a) 
2r [Ahsidj^hr) + Qh^ + Mos^^/il)] = ^2, (6.52b) 



where 

ai = TT{~g^{gfh2 - hg^) + {gfh^ - h^-gtYg^ - ~g^igfh^ - hi~gt) - {9fh2 - ^-9^)91}, (6.52c) 

and 

a2 - Tr{{T'~g^igfh2 ~ h^g^) + {gfh, ~ h^~gt)~g^ - g^{~gfhi - h,g^) - {gfh2 ~ h2~9^)gt)} ■ (6.52d) 

(6.52e) 

We note that although the boundary conditions of Kupriyanov and Lukichev are valid for large r, they also give the 
correct boundary conditions (that of continuity of the Green's functions and distribution functions) in the limit of 
r — > 0. For arbitrary transmission of a barrier with n channels, the boundary condition can be represented asi? 

- a - ^^S^nrr [gsl,gs2] fa ro\ 

gsAgsi = 7— V 2r„ . ^ . rr, (6.53) 

TT ^ A + in[gsigs2 + gs2gsi-2) 

where 7 is a constant factor, and T„ is the transmission of the nth channel. It can be seen that for a single channel 
with small transmission T, this equation reduces to the boundary condition of Kupriyanov and Lukichev. 
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7. PARAMETRIZATION OF THE QUASICLASSICAL GREEN'S FUNCTION 

The normahzation Ean. H6.11a|l permits a parametrization of the quasiclassical Green's functions that is very conve- 
nient for calculations. Ean. H6.11a|l is a matrix equation that is equivalent to the three equations H5.19|l for the s-wave 
component of the Green's function. To take into account the macroscopic phase of the superconductor, we note that a 
gauge transformation that transforms the vector and scalar potentials according to Eans. ()3.20|l and H3.21|l transforms 
the field operators according to the equations 

^ ^ i^e^^ (7.1a) 
i}+^^+e-^^ (7.1b) 

The Nambu-Gorkov Green's functions defined in Eans. H4.7|l are transformed accordingly. For example, two compo- 
nents of would transform according to 

[Gi2\ii ^ [Gf2"]iie-*(^("^i)-^(''^» (7.2a) 
[Gi2\i2 [Gf2"]i2e*^^^"^^+^^"^" (7.2b) 

Making the transformation (|2.25a|l to mixed coordinates, and taking the limit as r — > 0, we see that the off-diagonal 

components of the Nambu-Gorkov Green's functions are multiplied by a phase factor e'^^'^^'^ or e~'-^^^\ while the 
diagonal components remain unchanged. Consequently, A also transforms as 

A ^ Ae'x(^) (7.3) 

Keeping this in mind, we can express as 

~B. _ ( cosh6' sinh6'e'^\ cj a\ 

^ \- sinh 6ie-'x _ cosh 6 ) ^ > 

where 9 and x ^'^^ complex functions of the energy E and position R. This form satisfies the normalization condition 
gf- = t'^M' For completeness, we also give the expression for gf 

_^_/-coshr -sinh6l*e*x*\ 

~ i^sinh6l*e-*>^* cosh 61* / ^ ' 

We now put this into the Usadel equation for g^, Ean. (|6.22a|l . Keeping in mind that the matrix for A involves 
additional factors of e''^ and e~^^ due to the gauge transformation, the diagonal (1,1) component of this matrix 
equation is 

D sinh^ 9 d\x + D sinh 261 dgx ^§9 - 2i3( A) sinh 9 = 0, (7.6a) 
and the off-diagonal component (1,2) is 



Dd%9 - y sinh 261 (9^x)^ + 2Ei sinh 9 - 2i3fi( A) cosh 9 = 0, (7.6b) 



where we have used Ean. H7.6a|l to simplify Ean. (|7.6b|l . Defining a current js{E, R) by the equation 

js{E, R) = sinh^ 9{E, R)dj^x[E, R), (7.7) 

we can rewrite Ean. (|7.6a|l as 

Dd^jsiE, R) - 2i3?(A) sinh0 = 0. (7.8) 

js{E, R) is proportional to sinh^ 9, which is proportional to the square of the pair amplitude, and it is also proportional 
to the gauge-invariant gradient of the phase. Consequently, it is similar in form to the conventional definition of the 
supercurrent, and is called the spectral supercurrcnt. Indeed, Ean. l7.8|l is simply another way of writing Ean. l6.40|l . 
and Q and js are related by 

Q{E,R) = -^{j,[E,R)). (7.9) 
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As we noted before, dj^Q = if A is purely real, and from Ean. (|7.8() . it can be seen that dj^js{E, R) = also if A is 
real. 

Eauations H7.6|l form a set of coupled equations that can be solved in principle for 0{E, R) and x{E, R)- In the case 
of a negligible spectral supercurrent, the equations decouple, and one needs to solve only Ean. (|7.6b|) . In the limit 
of a bulk superconductor with a uniform real order parameter and no phase gradient, we recover the bulk value of 
the Green's function, Ean. H5.32|l . The differential equations must be supplemented by boundary conditions. From 
Ean. (|5.32|l . in a superconducting reservoir, we have 

cosh 6*50 = , (7.10) 
so that the value of 9 in the superconducting reservoir is given by 



? + ^ln^ ifi?<|A| 



^ 21^1^ ifi?>|A| 

The value of x in a superconducting reservoir is just the macroscopic phase of the superconductor. In a normal 
reservoir, 9 = 0. The value of x in a normal reservoir is meaningless, of course, and any choice that results in no 
phase gradient is valid. 

In terms of the 9 parametrization, the boundary conditions of Kupriyanov and Lukichev can be expressed as 

r sinh6'i(9^Xi) = sinh6'2 sin(x2 - Xi)j (7.12a) 

and 

r [9^6'i +isinh6li cosh 01 (9^X1 )] = cosh 6*1 sinh e'ae'^^^-'^i) - sinh6li cosh 6*2. (7.12b) 

Note that for r = 0, the boundary conditions reduce to xi = X2 and 9i = 02- In the absence of a supercurrent, 
{dj^Xi) = 0, so that the equations above simplify to 

Xi = X2, (7.13a) 

and 

r{dj^9i) = sinh(6l2 - 9i). (7.13b) 

Finally, we can write expressions for physical quantities in terms of 9 and x- These quantities can be written in 
terms of the Mij 

Mao = ^ [1 + cosh 6* cosh r - sinh6isinhr cosh(23(x))] , (7.14a) 
M33 = ^ [1 + cosh 61 cosh r + sinheisinhr cosh(23?(x))] , (7.14b) 
Mo3 = i sinh 61 sinhrsinh(25(x))), (7.14c) 
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and 



M30 = -isinh6'sinhr sinh(23(x)). (7.14d) 

If X is real, then the Afo3 — Aho = 0, and Mqo and A/33 simplify to 

Moo = cos2(3(6')) (7.15a) 

and 



M33 = cosh2(5R(6')) 



(7.15b) 
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These relations will be used in the next section when we discuss applying the Usadel equation to derive the transport 
properties of some simple device geometries. 

To conclude this section, we shall write an expression for the gap in terms 0. Replacing Ean. (j5.20|l for into 
Ean. (|5.13|l for the gap, we have 

A = No^ldE [gfh - h~g^]^2, (7.16) 

where performing the angular average in Ean. (|5.13|l gives the s-components of the Green's functions. Putting in h in 
the form of Ean. (|5.23|l . we obtain 

A = No^J dE [hL{~gf - ~gf) + hrigfr' - T^~gf)U (7.17) 

With gf and gf given by Eans. (|7.4|) and l|7.5|l . we obtain 

A^iVo^ j dE /iL(sinh6le^'^+sinhre'x*)-/iT(sinh6'e*X-sinh6'*e^'^')) . (7.18) 

As an example, consider the case of a bulk superconductor, where x = and Ht = 0. We then have 

A = A^o^ j dE /iL5R(sinh6i) 

= iVoA ^ dE tanh(S/2fcBT)5R (^^=A_ 

which is the usual self-consistent equation for the gap. 



(7.19) 



8. APPLICATIONS OF THE QUASICLASSICAL EQUATIONS TO PROXIMITY-COUPLED SYSTEMS 

We shall conclude our discussion of the quasiclassical theory by applying the equations that we have derived to 
some simple devices incorporating normal metals in close proximity with superconductors. Since the equations of 
motion in the diffusive limit are in general nonlinear, solving them usually requires numerical techniques, except in 
the limit of large resistances between the superconductor and the normal metal, where the Usadel equation can be 
linearized. We shall restrict ourselves to one-dimensional examples; these are the ones discussed most in the literature. 



A. Proximity-coupled wire 

We start with the simplest possible device, a one-dimensional normal-metal wire of length L connected on one end 
to a superconducting reservoir and at the other end to a normal metal reservoir. For definitiveness, let us take the 
superconducting reservoir to be at a; = and the normal-metal reservoir at x ^ L, and let us consider first the case 
of a perfect SN interface, so that the interface barrier resistance parameter r — 0. In this geometry, there can be 
no supercurrent, so that Q (or alternatively, js) is zero. Furthermore, we can take the phase x to be zero in the 
superconducting reservoir without loss of generality, and we note again that A = in the normal metal. Under these 
conditions, only Mqq and are non-zero in Eans. (|(j.34|l . which now read 



and 



MooidM ^ Ki, (8.1a) 



Mssidj^hr) = K2, (8.1b) 
where Ki and K2 are constants of integration. On integrating these equations from a; = to x = L, we obtain 

hL{x = L) - hL{x ^ 0) ^ Ki [ -^dx, (8.2a) 

Jo ^Woo 

hrix ^ L) - hrix ^ 0) ^ K2 [ —dx. (8.2b) 
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To calculate the conductance of the normal-metal wire in the linear approximation, we apply a small voltage V on 
the normal-metal reservoir, keeping the superconducting reservoir at = 0. If we consider the second equation, 
hxix = 0) = 0, and expanding hrix = i) in a Taylor's expansion to first order, we obtain from Ean. l)8.2b|l 



eV 



2fcBTcosh^ 



1 



M33{E,x) 



-dx 



.3) 



The electric current in the linear response regime can then be obtained from Ean. H().35|l 

H -1 



3 



2kBT 



dE- 



1 



cos\i\El2kBT) 



1 



M33{E,x) 



dx 



(8.4) 



There are two differences between this equation and the equivalent equation for a normal metal in the classical regime 
fEan. (|1.10a|l that we derived earlier. First, the equation above involves an integral over energy and position. One 
can define a energy and position dependent electrical diffusion coefficient 



D3{E,x) ^ DM33iE,x) 



.5) 



instead of the constant diffusion coefficient D for Ean Hl.IUa|) . Second, the current in Ean. (|8.4|) does not involve 
temperature differentials. Indeed, if one assumes that the superconducting reservoir is at a temperature T, and the 
normal reservoir at a temperature T + AT, and expand hx^x = L) to first order in AT, the terms involving AT cancel, 
so that there is no term proportional to AT. As a consequence, thermoelectric phenomena cannot be described using 
the quasiclassical equations, and an extension to the theory is required to take them into account iS 
From Ean. H8.4|l . one can define a spectral or energy dependent conductance of the wire 



G{E) = G 



N 



M33{E,x) 



dx 



where Gn is the normal state conductance of the wire. The total conductance is then 



G^ dE 



G{E) 



2kBTco&\i{E/2kBT) 



(8.6) 



■7) 
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FIG. 4: Spectral conductance G[E) of a one-dimensional wire of length L as a function of energy E, normalized to the 
correlation energy Ec- The barrier interface parameter r=0, and the gap is set to be A = IQQQEc- 
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Figure 4 shows the results of a numerical calculation of G{E)/Gm as a function of the normalized energy E/Ec- The 
normalization factor E^. — D/L^ is called the correlation energy or Thouless energy (from the theory of disordered 
metals, where it also occurs), ■^■^ and is dependent on the length L of the wire. At high energies, G{E) approaches 
its normal state value, as expected. As the energy is lowered, the conductance increases, as might be expected in 
a proximity-coupled normal metal. However, instead of continually increasing as the energy is reduced, it reaches 
a maximum of around 1.15 Gn at an energy of ~ 5Ec, and then decreases, reaching its normal state value at 
E = 0. This nonmonotonic behavior of the conductance is called reentrance, and has been observed in experiments 
by a number of groups. '^'^ It should be emphasized that the relevant energy scale where the maximum in conductance 
is observed is set not by the gap A of the superconductor, but by Ec, which itself depends inversely on the square 
of the length of the sample L. Hence in very long or macroscopic samples, the energy (and correspondingly, the 
temperature) at which the minimum would occur is far below the experimentally accessible range, and one regains 
the monotonic behavior expected from the simple Ginzburg-Landau theory of de Gennes.^* 

The temperature dependent conductance G{T) can be obtained from G{E) using Ean. H8.7|l : the result of this 
calculation, plotted in terms of the normalized resistance R{T)/Rn, is shown in Fig. 5. In obtaining this plot, we 
have used a value of A = 32Ec, corresponding to a weak-coupling transition temperature of Tc = 1.764A/fcB, typical 
parameters for Al films. We have also assumed that the gap is temperature dependent. Like G{E), R{T) is also 
non-monotonic, with a minimum at some intermediate temperature. We would expect the minimum in R{T) to be 
around T ~ bEc/ks, based on the behavior of G{E). However, the temperature dependence of the superconducting 
gap modifies this behavior, so that the minimum in resistance occurs at a somewhat higher temperature when the 
interface between the normal metal and the superconductor is perfect. Figure 5 shows additional curves corresponding 
to progressively increasing values of the interface barrier parameter r. Increasing the resistance of the NS interface 
decreases the leakage of superconducting correlations from the superconductor into the normal metal, and consequently 
results in a smaller increase in the conductance of the proximity-coupled normal metal. In addition, the temperature 
Tmin at which the minimum in resistance occurs is also shifted down. 
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FIG. 5: Temperature dependent resistance R{T) normalized to the normal state resistance Rn, as a function of the temperature 
T normalized to Ec, for different values of the interface resistance parameter r. The gap is set to be A = 32Ec. 



From Ean. Ht).38|l . the normalized density of states N{E) can be expressed in terms of 9 as 

N{E) = cosh {^{9)) cos (3(6*)) . (8.8) 

Figure 6 shows the density of states as a function of energy and position along the wire of length L. There is a 
proximity-induced decrease of N(E) near the superconducting reservoir. In fact, at the NS interface, there is a 
divergence in N{E) at the gap energy, and it goes to zero at £' = 0, just as one would expect for a superconductor. 
However, unlike a superconductor, it is not strictly zero for E < A, but still has a finite amplitude. As one moves 
away from the NS interface into the proximity-coupled normal wire, both the amplitude of this effective gap and the 
divergence are smoothly reduced, so that at the normal reservoir, one recovers the normal density of states. This 
position dependent variation of the density of states has been observed in experiments^ 
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FIG. 6: Density of state N{E, X) normalized to the normal state density of states A'o for a one dimensional wire of length L, 
as a function of E/Ec and position x along the wire. The superconducting reservoir is at a; = 1.0, and the normal reservoir is 
at a; = 0. The density of states is suppressed at low energies near the superconducting reservoir. 



In our analysis above of the proximity effect in a normal metal coupled to a superconductor, we have ignored the 
effects of electron decoherence on the proximity correction. Phase coherence is essential to observing the proximity 
effect; if the phase coherence length is less than the length of the sample L, a finite spatial cutoff of the proximity 
effect is introduced. Phenomenologically, this can be taken into account by saying that the length of the wire is now 
instead of L. Since is typically of the order of a few microns even at low temperatures, this sets the dimensions 
of the samples that are required to observe this mesoscopic proximity effect. 

While sets the upper cutoff for observing the proximity effect, a second relevant length scale for the problem 
can be obtained by considering the length at which Ec is equal to ksT. This length is 



hD 



(8.9) 



where Lt is called variously the thermal diffusion length or the Thouless length, again from the theory of disordered 
metals, where it also occurs. "^^ (We put in here explicitly the factor of h.) In fact, here it is simply the diffusive form 
of the superconducting coherence length in the normal metal, familar from the de Gennes/Ginzburg-Landau theory 
of the proximity effect)^ which in the clean limit is given by 



N 



(8.10) 



At low temperatures, when Lt is longer than L, the superconducting correlations induced in the normal metal extend 
throughout its length. At higher temperatures, they are restricted to a region of length Lt near the superconductor. 
Lt is also on the order of a few microns in typical metallic samples in accessible temperature regimes, so it also sets 
a limit on the dimensions of the samples in which one can see a proximity effect. 

To calculate the thermal conductance of the wire, we proceed from Ean. (|8.2a|l . Wc now consider a small temperature 
differential AT applied across the wire. Expanding in a first-order Taylor's expansion, as we did for Ht, we obtain 



We then obtain from Ean. Hfci.36|) 



SAT 



Mqo{E,x) 



-dx. 



(8.11) 



Jth = — -- / aE- 



2kBT^ J cosh' iE/2kBT) 



1 



Moo{E,x) 



-dx 



.12) 
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As with the electrical conductance, we can define a thermal difiiision coefiicient 

DoiE,x) = DMooiE,x), (8.13) 

and a spectral thermal conductance 



Gth{E) = GthN 



- 1 



-dx 



(8.14) 



/o MooiE,x) 

where GthN is related to the normal state electrical conductance by Eon. 1)1. 17(1 

GthN = Gat— — ^-T. (8.15) 
3 

Finally, the thermal conductance itself is given by an integral over energy 

= ^^(w/^'^cosh^l'iJ/IST)- ^'-''^ 

Of course, as noted by Andreev^^ the thermal conductance of a normal metal wire sandwiched between a normal-metal 
reservoir on one end and a superconducting reservoir on the other end must vanish, since the superconductor acts as 
a thermal insulator, so that no thermal current can flow through the device as a whole. However, one may consider a 
normal-metal wire with the superconducting reservoir connected off to one side, so that it does not block the flow of 
thermal current through the proximity-coupled wire. One may then consider the thermal conductance of the normal 
metal wire itself. Figure 7 shows the thermal conductance of this geometry, as a function of temperature, for different 
transmissivities of the interface barrier. The thermal conductance shows a monotonic decrease as T is lowered below 
Tc, although there are no distinct features at any particular temperature, unlike for the electrical conductance. In a 
superconductor, the exponential decrease in the thermal conductivity is associated with the opening of the gap in the 
quasiparticle density of states, since it is the quasiparticles that carry the thermal current. Noting the decrease in the 
density of states in the proximity-coupled normal metal wire, shown in Fig. 6, it is not surprising that this system 
will also show a decrease in the thermal conductance. The thermal conductance of the wire is strongly dependent 
on the transmission of the NS interface, characterized by the parameter r, and approaches the normal state thermal 
conductance as r increases. 




FIG. 7: Thermal conductance Gth of a normal wire connected to a superconducting reservoir on one end, and a normal metal 
reservoir on the other, as a function of the normalized temperature T, for a number of different values of the interface barrier 
parameter r. The gap is set to be A = 32Ec, corresponding to a transition temperature of T = 18.14_Ec. 



We note here again that, in our current approximation, a small voltage drop across the S/n-wire/N device will 
not result in a contribute to a thermal current through the system, since any terms proportional to voltage in the 



33 



expansion of will cancel. This is the converse of the case for the electrical current, where a small temperature drop 
did not contribute to the electrical current, emphasizing again that the conventional quasiclassical approximation 
cannot take into account thermoelectric effects. 



B. Superconductor- metal bilayer 

Instead of a normal reservoir on one side of the wire, if we consider a wire of length L connected only on one 
end to a superconducting reservoir (with the other end open), then a true gap in the density of states opens up 
in the proximity-coupled normal metal. The magnitude of the gap is related to Ec] hence, one can consider the 
proximity-coupled normal-metal in this case as a superconductor with a gap of Ec. 

If the superconductor is not a reservoir, but a thin layer itself, then one will suppress superconductivity in the 
superconducting layer due to the proximity of the normal metal, an inverse proximity effect. The suppression of 
superconductivity is expected to reduce the gap in the superconductor. It is an interesting exercise to calculate the 
transition temperature of the bilayer in the quasiclassical approximation. We shall loosely follow here the treatment 
given by Martinis et alM- and Gueroni^ 

Let the thickness of the superconductor be ts, and the thickness of the normal metal t^- We take the origin, a; = 0, 
at the NS interface, and the superconductor extends from x = —ts to a; = 0, and the normal metal from cc = to 
X = tfq. Near the transition, the order parameter in the superconductor is small, so we may the small 9 limit of 
Ean. (|7.6b|) . Since the phase is not important in this problem, we take x = 0. The resulting equation is 

Ddl0 + 2Ei9-2iA = O, (8.17) 

where we have assumed that the gauge is chosen so that A is real. Let us assume that 9 at x ~ in the superconductor 
is 9qs, and that variations of 9 about this mean value are small. Under these conditions, we can also assume that the 
gap A in the superconductor is uniform. We can then expand 9s to second order in x 

9s = 9os + ax + bx^. (8.18) 

Now, at X = —ts, (and also at x = t^), we have a vacuum interface, where dx9 = 0. Hence a = 2hts, and from the 
differential equation (|8.17|) taken at ec = 0, & = {i/D){A — E9qs), so that 

Os - 0ns + ^(A - E9os){2tsx + x^). (8.19) 

Similarly 

i 

D 

From the boundary condition Ean. lj7.13b|) . we have the two equations 

2irts 



N = (^0N + —E9oN{2tNX- x^). (8.20) 



D 

Solving this pair of equations for 9os, we have 



(A - E9os) = 0ON - Oos, (8.21a) 
E9qn — Oqn ~ 9qs- (8.21b) 



D 

2irt]^ 



n _ A EHs^ts + In) + ^rH\tl - 2irDtl 



Putting this into Ean. (|7.19|l for the gap, with T = Tc, and with the approximation that sinh(0) ~ 9, we obtain 



1-iVo^ / ^ tanh(£;/2fcsrc 



^ tN{ts + tN) 



{ts + tNy + i'iryD^)titj, 



.23) 



The first term in the square brackets gives the bare transition temperature Tco of the superconducting film, and the 
second term corresponds to the corrections associated with the inverse proximity effect. For a perfect interface, with 
r = 0, the suppression of is directly proportional to the normal fraction of the bilayer, and T^, ^ as tj^ increases. 
For a highly resisitive barrier (r — > oo), the second term in the square brackets goes to zero, so that there is little 
effect of the normal metal film on Tc of the superconducting film, as expected. 
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C. The SNS junction and Andreev interferometers 



As our final example of the application of the quasiclassical equations of superconductivity, we consider the case 
of a dirty SNS junction. The model we consider is a normal metal wire of length L sandwiched between two su- 
perconducting reservoirs. As is well known, the application of a phase difference between the two superconducting 
reservoirs will result in the flow of a supercurrent through the normal wire. The phase difference can be applied, for 
example, by connecting one of the superconducting reservoirs to the other, thereby forming a loop with two different 
arms, one superconducting and one normal. This configuration is commonly called an Andreev interferometer. The 
phase between the two superconducting reservoirs can be varied by applying an Aharonov-Bohm type magnetic flux 
through the area of the loop; in this respect, we put together all contributions in the gauge-invariant phase x- Due to 
the single- valued nature of the wave functions, a phase factor of 27r<l'/<I'o is picked up in going aroung the loop, where 
<i> is the magnetic flux threading the Andreev interferometer, and $o = h/2e is the superconducting flux quantum. 
In a superconductor, the supercurrent Is that is generated is directly proportional to the phase gradient; if Ig is 
small compared to the critical current Jc, the phase dropped across the superconductor will be small. Since Ic of 
the superconducting part of the Andreev interferometer is so much greater than the critical current of the proximity 
coupled normal-metal wire, most of the phase change will occur across the length L of the normal metal wire. 

This fact allows us to map the Andreev interferometer that is coupled with an Aharonov-Bohm fiux <I> to a SNS 
system with a phase difference (j) = 27r$/$o across it. In terms of our model, we consider the superconductors 
to be reservoir; this means applying a boundary condition for the gauge-invariant phase x at the superconducting 
reservoirs. For our purposes, we apply this boundary condition anti-symmetrically, with a phase xl = — 7r$/$o at 
the superconducting reservoir at x = 0, and XR = 7r$/$o at the superconducting reservoir at x = L. We then must 
solve Bans. 17. 6() in the normal-metal wire for 6 and Xj with A = 0, subject to the boundary condition for x noted 
above, and the boundary condition 9 = Oso (where ^50 is given by Ean. H7.11|l ) at both x — Q and x ~ L. In general, 
both 9 and x are complex functions of x and £', and the solution of Eans. (|7.()|l must be done numerically. 
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FIG. 8: (a) Spectral supercurrent Q in a normal wire between two superconducting reservoirs, as a function of the normalized 
energy E/Ec, for different values of the phase difference tj) between the superconducting reservoirs, (b) Geometry of an Andreev 
interferometer, essentially a cross with one arm connected to superconducting reservoirs, and the other arm connected to normal 
reservoirs. 

Following Yipj2ii we consider first the supercurrent Q, given by Ean. (|7.9|l . Some insight into the contributions to 
js can be gained by looking again at the case of a bulk superconductor. From Ean. Ht).46|) . the major contribution 
to js comes from energies near the gap. For a long proximity wire with Ec < A, however, the major contribution 
comes from energies of order Ec- Figure 8(a) shows a plot of Q as a function of energy for various values of the 
phase difference (p between the two superconducting reservoirs, for the case of a perfectly transparent interface, and 
A = 32Ec. Of course, for zero phase difference, the supercurrent vanishes. As <j> is increased from zero, there is a peak 
in Q{E) at £^ ~ QEc- This peak moves down in energy as (p increases. At larger values of E, Q becomes negative. 
For shorter wires, this region of negative Q is less prominent. The total supercurrent is given by the second term in 
Eqn.(|n23 

Js = eNoD J dE Q{E)hL{E), (8.24) 
and therefore depends also on the distribution of quasiparticles. Any change in this distribution will affect the 
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supercurrent. For example, the distribution can be changed by increasing the temperature, which has the result of 
decreasing the supercurrent. As we demonstrated in the introduction, a non-equilibrium quasiparticle distribution 
can also be generated by injecting a normal current into the proximity wire in the SNS geometry by attaching 
two additional leads to the center of the normal wire, forming a normal cross, with two of the wires attached to 
superconducting reservoirs, and the two other wires attached to normal reservoirs, as shown in Fig. 8(b). The current 
in the SNS junction is then a function of the current injected between the normal reservoirs, and the supercurrent 
can even change sign depending on magnitude of the injected normal current^SiiJi This effect has been observed in 
recent experiments4i 

Due to long-range phase coherence, the Green's function in the arms of the cross attached to the normal reservoirs 
will also depend on the phase difference between the two superconducting reservoirs in the structure shown in Fig. 8(b). 
Consequently, the electrical conductance measured between the two normal reservoirs will also be a periodic function 
of the phase difference between the two superconducting reservoirs. Experimentally, the conductance of such Andreev 
interferometers have been found to oscillate periodically with an applied external flux, with a fundamental period of 
4>o — /i/2e42iiii4 Similar oscillations are also expected in the thermal conductance as well. Periodic oscillations are 
also observed in the thermopower of Andreev interferometers,'*^ although these thermopower oscillations cannot be 
described within the framework of the current quasiclassical theory. 



9. SUMMARY 



The quasiclassical theory of superconductivity has proved to be a powerful tool for the quantitative description of 
long-range phase coherent phenomena in diffusive proximity coupled systems. As we have shown, the linear electrical 
and thermal conductance of complicated devices incorporating normal and superconducting elements can be calculated 
in principle, although the solutions frequently involve numerical techniques. Extension to the nonlinear regime, with 
finite voltages across the normal reservoirs, is also conceptually straightforward, although numerically challenging. 

Application of finite voltages to the superconducting elements is trickier, as it involves time dependent evolution of 
the phase, and is only beginning to be examined theoretically. Finally, the quasiclassical theory for diffusive systems, 
in its present form, does not deal at all with thermoelectric phenomena. Extensions to incorporate thermoelectric 
effects in the theoretical framework have been attempted^^ but still require further work to be complete. 
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We set ft = 1. 

From the right, we operate on the coordinate 2, hence we use the operator Gq2 ■ operated on the right is equivalent to 
Gq2* operated on the left. 

A defined this way differs from the conventional definition by a factor of i. Note that the signs associated with A in the 
matrix above are different from the conventional definition of Gorkov's equation, because of the factor in the definition 
of the Green's functions. 

^" Since we are dealing with the static limit in all that follows, we again use the symbol T to refer to the temperature 
One can also express gf equivalently in terms of sines and cosines. 
See, for example, F. Wilhelm, PhD Thesis, Universitat Karlsruhe, 2000 



